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^ '■ ABSTRACT 

Oh' We present deep ground based imaging of the environments of five QSOs that contain 

q , sub-Damped Lyman-a systems at z<l with the SOAR telescope and SOI camera. 

5h ■ We detect a clear surplus of galaxies in these small fields, supporting the assumption 

ryj ' that we are detecting the galaxies responsible for the absorption systems. Assuming 

& . these galaxies arc at the redshift of the absorption line systems, we detect luminous 

L>L* galaxies for four of the five fields within 10"of the QSO. In contrast to pre- 
vious imaging surveys of DLA systems at these redshifts, which indicate a range of 
. morphological types and luminosities for the host galaxies of the systems, the galax- 

(*T) ' ies we detect in these sub-DLA fields appear to be luminous (L>L*). In the case of 

£SJ , the absorber towards Q1009— 0026 at z a f, s =0.8866 we have spectroscopic confirmation 

CNJ ■ that the candidate galaxy is at the redshift of the absorber, at an impact parameter 

' of ~35 kpc with a luminosity of 3 < L/L* < 8 depending on the magnitude of the 

A^-correction. These observations are in concordance with the view that sub-DLAs 
■ may be more representative of massive galaxies than DLA systems. The environments 

\ of the absorbers span a range of types, from the inner disk of a galaxy, the periphery 

of a luminous galaxy, and the outskirts of interacting galaxies. The large impact pa- 
rameters to some of the candidate galaxies suggest that galactic outflows or tidal tails 
arc likely responsible for the material seen in absorption. Wc find a weak correlation 
between Nh i and the impact parameter at the 2er level, which may be expected from 
^ ' the heterogeneous population of galaxies hosting the absorption line systems and ran- 

dom orientation angles. In addition, we detect a possible gravitationally lensed image 
of the BL-Lac object Q0826-2230. 
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1 INTRODUCTION 

One of the principal manners in which to study high redshift 
galaxies is through the use of Quasar (QSO) absorption line 
systems. These absorption line systems with the largest Nh i 
as measured by the Lyman-a A 1215 line are classified into 
two groups, the sub-Damped Lyman-a systems (sub-DLAs: 
19.0 < log Nh i < 20.3 cm -2 ) and Damped Lyman-a sys- 
tems (DLAs: log N H i > 20.3 cm -2 ). Although it is rela- 
tively easy to obtain high S/N spectra of QSOs with the 
large telescopes and sensitive high resolution spectrographs 
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available today, much of the difficulty in interpreting the 
data from spectroscopic observations of QSO absorbers is 
that the environment in which the absorbing gas is resid- 
ing is almost always unknown. There have been many stud- 
ies aimed at imagin g of DLAs (e.g.. iLe Brun et al.lll997l: 



Kulkarni et ail I2000I. l200ll: IWarren et al, . 
20031 : Christensen et al, 1 12004 . l2005l ; IChen et al.ll2003l . l200 



200lFTRao et al 



However, many of these studies have resulted in non- 
detections, or detections of dwarf or low surface brightness 
galaxies. Higher redshift systems (z > 1) are increasingly 
difficult to detect, with only ~10 percent of the host galax- 
ies of DLAs identified. 

The non-detection of luminous galaxies in most DLAs 
searched for, together with the low metallicities and weak 
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metallicity evolution determined via the absorption lines 
of undeple ted elements such as S o r Zn in the DLA sys- 
tems (e.g. IProchaska fc Wolfe! 120021 : iKulkarni et"aD 120051 : 



iMeiring et all l2006t ) suggests that while these objects 



do contain the m ajority of neutral gas in the Universe 
(jWolfe et al.lll995h . they may be tracing mainly star for- 
mation in lower mass galaxies. Chemical evolution models 
of dwarf galaxies predict low metallicity that slowly evolves 
until more recent times (|Pagel fc Tautvaisienelll998h . 

We have recently completed a survey of 31 sub- 
DLAs at z < 1.5 using the Magellan Inamori Kyocera 
Echelle ( MIKE) and Ultraviolet Visual Ech el le Spectrograph 
(UVES) (iMeiring et al.ll2007l , |200S] , 12009 A| |b1; IKulkarni etall 
120071 ; IPeroux et all l2006al lbl. (200a ). At odds with the typi- 
cally low metallicities seen in DLA systems, several of the 
sub-DLAs in the sample showed near solar or above solar 
metallicity (i.e. [Zn/H]>0.0). In fact, a trend towards in- 
creasing metallic ity with decreasing Nh i is seen over thr ee 
decades in N H 1 (|Meiring et al.ll2009Al ; Ikhare et al.ll2007l ). 

Several possibilities exist to explain such a trend: 1) 
there is a dust obscuration bias that limits the num- 
ber of observable metal rich DLAs to b e observed as the 
amount of dust increase s with metallicity (|Boisse et al.ll 199*81 ; 
IVladilo fc Perouxll2005l ); 2) the low H I column densities in 
the sub-DLA regime are a result of conversion of neutral hy- 
drogen into molecular gas, which is forming stars and pro- 
ducing the metals; 3) the lower H I column densities are 
a result of increased ionisation from star forming regions 
that produce metals that enrich the surrounding areas; 4) 
sub-DLAs are intrinsically more likely to aris e in massive 
galax ies which are typically more metal rich (|Khare et al.l 
120071 ). In all likelihood, it is possible that all of these effects 
are at work together to produce the observed trend. 

Few su b -DLA s have been direct ly imaged to date 
llChen et alj 120051 ; IStraka et all l20ld ; IChristensen et all 
l2005h . With this in mind, we have started a program to 
conduct follow up imaging of sub-DLAs at z < 1.0 to 
determine the nature of these systems. Here, we present 
deep SOAR imaging of the five QSO fields Q0826-2230, 
SDSS J1009-0026, SDSS J1228+1018, SDSS J1323-0021, 
and SDSS J 1436-0051. The structure of this paper is as fol- 
lows: the observations are discussed in § 2, notes on the 
individual fields are given in § 3, results of the observations 
are given in § 4, and a discussion is given in § 5. Spectra 
of the QSOs showing the absorpt i on lines can be found in 
IMeiring et all |2007i , |2008| , 1200961 ): IPeroux et all |2006al ). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

These data were obtained with the 4.1m Southern Observa- 
tories Astrophysics Research (SOAR) telescope located at 
Cerro Pachon, and the Soar Optical Imager (SOI) camera 
during February and March 2009. Roughly 1.5 nights were 
lost due to weather during the 4 night run. The SOI camera 
has a pixel scale of ~0.154 arcsec pix - with 2x2 binning, 
which was used throughout the observations. The two 2kx4k 
detectors are arranged in a mosaic configuration with a ~8 
arcsec gap between the detectors. The total field of view of 
this instrument is 5.25' x 5.25'. 



Each field was observed in some combinatio n of the 
Sloan Digital Sky Survey (SDSS) l|York et al.ll200ol ) g, r, i, z 
filter set. The targets in this investigation were selected 
mainly due to visibility at the SOAR telescope and the red- 
shift of the absorber, favoring the lower redshift systems in 
our prior spectroscopic sample to increase the chances of 
detections. Exposure times in each band and coordinates of 
the background QSO are given in Table [1] Each exposure 
was dithered by ~10 arcsec to fill the gap between the de- 
tectors and minimize the effects of bad pixels. Individual 
exposure times ranged from 300 to 900 seconds, depending 
on the magnitude of the background QSO. Series of twilight 
sky flats were taken at the beginning of each night for each 
filter and used for the flat-field images. 

We have reduced these data 

with the theli pipeline available at 
http : / / www . a stro . uni-bonn . de/ ~mischa/ theli .html 
jErben et all 120051 ). theli is a multipurpose pipeline 
configurable for many mosaic imagers. Data reduction 
proceeds with typical steps for imaging data. The mosaic 
frames are first split, then overscan corrected and trimmed. 
Next, the frames are bias corrected. A flat field image is 
made from the twilight exposures taken with the same filter 
used for the science frames. The data are then processed 
and a superflat is created from the smoothed data and a 
fringe model is also created for fields observed in the i and 
z bands. The fringe model is then scaled and subtracted 
from the calibrated data. 



We used the SCAMP package for astrometric calibration 
and SWARP for combining the individual frames. Both of 
these packages can be found at http://terapix.iap.fr/. 
The Sloan Digital Sky Survey (SDSS) astrometric catalog 
was used for all fields except for Q0826— 2230, for which we 
used the 2 Micron All Sky Survey (2MASS) for astrometric 
references. The fields were photometrically calibrated from 
the SDSS catalogues for the fields with multiple sources 
with m < 20.0 detected in both the SDSS images and these 
data. The photometric zero point for the field of Q0826- 
2230, which is outside of the SDSS footprint was calibrated 
by determining the slope of the best fit line between the 
photometric zero points and airmasses of these fields as well 
as four others not published here. SEXTRACTOR version 2.5.0 
was used for object iden tification and final catalog creation 
jBertin fc Arnold 1 19961 ). A 7x7 convolution mask of a 
gaussian PSF with FWHM of 4 pixels was used to filter the 
images during object detection. Only objects greater than 
3a above the background noise were included in the final 
catalogs, with a minimum area of ten pixels. Magnitudes 
given in the text are isophotal, and centroid positions 
are windowed. The magnitudes have been corrected for 
Galactic extinction using the online calculator available at 
http : //www. ipac . caltech.edu/forms/ calculator .html. 
All magnitudes given here are in the AB system. 

Three band colour images we re also made of the f rames 
using the prescription detailed in iLupton et all |2004) . The 
images were aligned in pixel space, and a scaling was applied 
to each channel that attempted to match the colours of the 
objects seen in the SDSS images. 
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Table 1. Table of observations. 



Q0826-2230 






Zem > 0.911, a = 08 


: 26 : 01.5, S = 


-22 : 30 : 


26.2 




9 


r 


i 


z 


Integration Time (s) 


600 


3300 


3000 


5400 


FWHM (") 


1.0 


1.0 


0.7 


1.25 


Limiting Mag. 


24.5 


25.0 


25.0 


23.5 


SDSS J1009- 


-0026 






z em = 1.241, a = 10 : 


09 : 30.40, 5 = 


-00 : 26 


: 19.1 




9 


r 


i 


z 


Integration Time (s) 


6900 


3600 


5400 


2700 


FWHM " 


0.9 


0.8 


1.1 


1.1 


Lim. Magnitude 


26.0 


25.5 


25.0 


23.5 


SDSS J1228+1018 






z em = 2.306, a = 12: 


28 : 36.80, 5 = 


+10 : 18 


: 41.7 




g 


r 


i 


z 


Integration Time (s) 




3600 


3600 


8100 


FWHM " 




0.8 


0.8 


1.5 


Lim. Magnitude 




25.0 


24.5 


23.0 


SDSS J1323- 


-0021 






Zem = 1.388, a = 13 : 


23 : 23.78, 5 = 


-00 : 21 


: 55.3 




g 


r 


i 


z 


Integration Time (s) 


2100 


3600 


4500 




FWHM " 


0.7 


0.8 


0.7 




Lim. Magnitude 


25.5 


25.5 


25.0 




SDSS J1436- 


-0051 






z em = 1.275, a = 14 : 


36 : 45.03, 5 = 


-00 : 51 


: 50.6 




g 


r 


i 


z 


Integration Time (s) 


7200 


3600 


4800 


3600 


FWHM " 


0.8 


0.8 


0.6 


1.4 


Lim. Magnitude 


26.5 


25.5 


25.0 


24.0 



2.2 PSF Subtraction 

Galaxies at small radii to the background QSO may become 
visible after PSF subtraction. With the 5' x 5'field of view of 
the instrument, several suitable point sources are available 
in each image for reconstructing the point spread function 
(PSF) of the instrument. The PSF did not vary significantly 
over the field of view of the instrument. Factors such as 
slighly different astrometric solutions and distortions deter- 
mined for the individual frames as well as tracking jitters 
due to wind can lead to non-ideal PSFs in the final com- 
bined frames. We have used the daophot package in IRAF 
to create model PSFs using multiple stars in the images. 
The Moffat profile was used with /3 = 2.5, which provided 
the best fitting model of the PSF compared to other options. 
A 30 pixel wide aperture was used in all cases. 

This model PSF was then used to inter- 
actively subtract the QSO in the images using 
the IDP3 package in IDL. IDP3 is available at 
http://mips.as.arizona.edu/MIPS/IDP3/ for down- 
load. The PSF is interactively shifted and scaled to 
optimally subtract the object of interest. 



In the case of Q0826— 2230, the peak brightness of the 
QSO itself in the r and i bands was well above the level 
where the CCDs in the mosaic preform linearly. As such, 
the PSF needed to be slightly over-subtracted for a proper 
fit in these images. Accurate PSF subtraction however was 
not a main goal of these observations, rather deep imaging in 
multiple bands to allow for galaxy colours to be determined 
was deemed a more important intent of these observations. 

3 COMMENTS ON THE INDIVIDUAL FIELDS 

3.1 Q0826 2230: z em >0.9110 2 a6s =0.9110 

Q0826— 2230 is a BL-Lac object with no known redshift 
due to the lack of any emission lines, although the absorp- 
tion system at z a h s =0. 9110 does giv e a lower limit to th e 
redshift of this object (|Falomd [l990l : iMeiring et ai1l2007l ). 
The sub-DLA in t he spectrum of this object has log Nh i 
= 19.04+0.04 (e.g. iRao. Turnshek fc Nestorl l2006t ). and a 
high metallicity of [Zn/H]=+0.68±0.08 with a kinematical 
width (defined as the width of the profile where the inner 
90 per cent of the absorption occurs) of Ati9o =278 km s -1 
based on our earlier high resolution spectra |Meiring et al.l 
120071 ). 

This field was observed with the SDSS g, r, i, z filters. 
Exposure times for each filter are given in Table [1] The 
final calibrated frames are shown in the top row of Figure [I] 
and the colour image composed of the g, r, and i images is 
shown in Figure [2] The seeing in the final combined images 
is ~1.0, 1.0, 0.7, and 1.25 arcsec in the g,r,i, and z filters 
respectively. 

Several objects are detected within 10 arcsec of the 
QSO, corresponding to ~78 kpc at z a()s =0.9110. As this ob- 
jects is near the plane of the Milky Way disk (I =243.99, 
b =+8.93), the stellar density is high in the region. An ob- 
ject 1.4 arcsec to the west of the QSO is seen most clearly 
in the i band image. A Point Spread Function (PSF) model 
of ten stars near the QSO was constructed by fitting 2-d 
Moffat profiles to the individual stars and an average PSF 
was created that was used to subtract the QSO from the i 
band image. The PSF subtracted i band image is shown in 
last column of Figure [T] 

Limiting magnitudes of ~ 24.5, 25, 25 and 23.5 were 
reached in the g, r, i, z filters in the exposure times given 
in Table Q] Objects labelled four and five in Figure [T] are 
the only possible objects that could be galaxies at z~0.9, 
objects three and six are certainly stars. Assuming object 
3 is at the redshift of the absorber, L>10L* making it un- 
likely to be the host galaxy of the absorption system due 
to the extre me luminosity and the paucity fo such galax- 
ies, hyperz (|Bolzonella. Miralles fc Pellal2000l ) was used to 
determine photometric redshifts for objects four and five of 
Zphot = 0.66 ± 0.10, 0.62 ± 0.14 respectively. Magnitudes for 
the objects are given in Table [2] 

Object labelled number two in Figure [TJ located ~1.5 
arcsec west of the brighter object is more clearly visible in 
the PSF subtracted image shown in right column of Fig- 
ure [T] There are several suggestions that this QSO may be 
gravitationally lensed. The colours of this object are quite 
similar to the main object, the QSO has a large luminosity 
(Mi < —28.2) indicative of magnification, and the side-by- 
side lensing pattern is the most commonly found lensing 
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patter n for QSOs (|Kochanek, Schneider PH 120041 ). iFalomol 
l|l990l ) also noted that this QSO is analogous to other lensed 
objects. Hubble Space Telescope imaging would likely reveal 
if this object is indeed gravitationally lensed. 



3.2 Q1009 0026: z em — 1.241 z abs =0. 8426, 0.8866 

There are two sub-DLAs in the spectrum of this QSO. 
One system has an absorption redshift of ,z a j, s =0.8426 with 
a low metallicity ([Zn/H]< -0.98 and [Fe/ H]=-1.28 and 



a kin ematical width of Avgp=36 km s iMe iring et al 



(120071) 1 and log N H i =20.20±0.06 (|Rao. Turnshek fc Nestor 
120061 '). and a second system at z ai)s =0.8866 has a high 
metallicity ([Zn/H]=+0.25±0.06 and [Fe/H]=-0.58 with 



A^9Q=94 kms~MMeiring et al l \200H) ) and log N H i =19.48 
(|Rao. Turnshek fc Nestorll2006h . 

This field was observed in the griz filters. Exposure 
times are given in Table [T] The final combined frames are 
shown in the second row of Figure [T] and the colour com- 
bined image composed of the gri frames is shown in Figure 
[5] The seeing in the final combined frames was ~0.9, 0.8, 1.1, 
and 1.1 arcsec in the g,r,i, and z filters respectively. Lim- 
iting magnitudes of ~ 26, 25.5, 25, and 23.5 were reached in 
the g, r, i, z filters. 

One object is detected in the images within 10 arcsec of 
the QSO, 4.6 arcsec to the west of the QSO (corresponding 
to a comoving impact parameter of ~35 kpc at z=0.8426 or 
~36 kpc at z=0.8866). This object has been confirmed to 
be at z = 0.8866 via integral field unit (IFU) observations 
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|Peroux et alj20ich . The hyperz software package was used 
to determine a photometric redshift of z p hot = 0.79 ± 0.08. 

We note that photometric redshifts for star form- 
ing galaxies are difficult due to the flatness of their 
spectrum, and four band photometry is a minimum re- 
quirement for determining photometric redshifts. Addi- 
tional photometric points in the Infra-red (IR) and near- 
UV (NUV) help to tightly constrain photometric redshifts 
(|Margonier. Connolly fc Tvsonl [20051 ) . No candidate galax- 
ies for the z a ;, s =0.8426 system are observed in the field. The 
limiting i band magnitude of ~ 25 corresponds to a limiting 
luminosity of this galaxy of L<0.1L* in the observed frame 
i band. No object at z ~ 0.8426 was observed in th e IFU 
observations of this QSO either (|Peroux et al.| [2 010b . The 
only other possible source in this field is a faint object south 
of the QSO which is minimally visible in the r and i frames, 
as well as the colour combined frame in Figure [2] The object 
was not detected at > 2a with SEXTRACTOR however, and 
was not included in the final catalogue. A second possibility 
is that the z — 0.8426 galaxy could be at a small impact pa- 
rameter to the QSO and not have significant star formation. 
The previous IFU observations would not likely detect such 
a galaxy, and higher quality PSF subtraction or space based 
imaging would be necessary to rule out such a scenario. 



image composed of the gri frames is shown in Figure [2] The 
seeing in the final combined frames was ~0.7, 0.8, and 0.7 " 
in the g, r, and i filters respectively. Limiting magnitudes of 
~ 25.5, 25.5, and 25 were reached in the g, r, i filters. 

Three objects are detected in the field, with one ob- 
ject detected in the i band image to the northeast of the 
QSO. After PSF subtraction of the QSO, the nearby object 
is more visible (labelled 2 in Figure [TJ and is in the same 
position as the object t hat was detected i n the K band IR 
imaging of this QSO in iHewett fe Wild I (120071 ) and adap- 
tive optics imaging of lChun et al. | (|2010h . The quality of the 
PSF subtraction of the QSO is not adequate to obtain a 
magnitude of this object in the i band image, and we note 
only its position in Table [2] 

It is inter esting to note that n o H-a emission is detected 
in this galaxy |Peroux et al.|[20fol ). The lack of H-a emission 
and the non-detections in the g and r images of this field 
(even though the seeing and depths were similar) indicates 
that this galaxy is of moderately early type, with few young 
stars producing emission in the near UV and u band at the 
redshift of the absorber. Indeed, IPeroux et all ( |2010h deter- 
mine a SFR of <0.1Mq yr _1 for this object based on the 
non-detection of H-a. Magnitudes for the remaining objects 
are also given in Tabled 



3.3 Q1228+1018: a em =2.306 z a6s =0.9376 

A sub-DLA system is observed in the spectrum of 
this QSO at z ai , 3 =0.9376 wi th log N H i =19.41±0.04 
jRao. Turnshek fc Nestol l200fj ). Although the Zn II A A 
2026,2062 lines were not detected in our earlier spectra, the 
system does have a high metallicity based on the depleted 
element Fe of [Fe/H] =— 0.31 ± 0.02 and a kinematical width 
Av 90 = l 16 km s" 1 jMeiring et alj|200gft . 

This field was observed in the riz filters with the ex- 
posure times given in Table [TJ The final combined frames 
are shown in the third row of Figure [T] and the colour com- 
bined image composed of the riz frames is shown in Figure 
The seeing in the final combined frames was ~0.8", 0.8", 
and 1.5" in the r,i, and z filters respectively. Limiting mag- 
nitudes of ~ 25.0, 24.5, and 23 were reached in the r, i, z 
filters. 

Two galaxies are seen northwest of the QSO at pro- 
jected impact parameters of 38 and 66 kpc with i band mag- 
nitudes of 22.08 and 22.70 respectively. Two other objects 
are also seen in the r band image of this field, but they were 
not detected above a 2a significance with SEXTRACTOR and 
are not included in the final object catalogues. 



3.4 Q1323 0021: z em = 1.388 z ai)s =0.7160 

A metal rich sub-DLA system is observed in the spectrum 
of this QSO at z ai; s =0.7160 with log N H i =20.21±0.20 
jPeroux et alj|2006ah . From high resolution VLT spectra of 
this object, metallicities based on Zn of [Zn/H]=+0.61±0.20 
and based on Fe of [Fe/H] =— 0.51±0.20 and a kinematical 
width Ai>9q=91 km s -1 h ave been determined (|Khare et al.l 
120041 ; IPeroux et af1l2006al ) . 

This field was observed in the gri filters with the expo- 
sure times given in Table [1] The final combined frames are 
shown in the fourth row of Figure[TJ and the colour combined 



3.5 Q1436 0051: z em = 1.275 ^^=0.7377,0.9281 

Two systems are observed in the spectrum of this QSO 
at Zabs =0.7377 and 0.9281 with log N H i =20.08±0.11 
and log Nh i <18.8 respectively. Both systems have high 
metallicities, with the lower redshift system at z a bs =0.7377 
having [Zn/H]=-0.05±0.12 and [Fe/H]=-0.61±0.11 with 
At>9o=71 km s -1 , while the system at z a (, 3 =0.9281 is 
highly enriched with [Zn/H]>+0.86 and [Fe/H]> -0.07 and 
Augo=62 km s _1 . Both systems show Ca II AA 3934 absorp- 
tion lines, with [Ca/Fe] = -0.98±0.03 and [Ca/Fe] = -0.79± 
0.04 fo r the Zai, 3 =0.7377 a nd z a (, s =0.9281 systems respec- 
tively jMeiring et alj|2008h . 

The field of this QSO was observed in the griz filters 
with the exposure times given in Table[TJ The final combined 
frames are shown in fifth row of Figure [1] and the colour 
combined image composed of the gri frames is shown in 
Figure [2] The seeing in the final combined frames was ~0.8, 
0.8, 0.6 and 1.4 arcsec in the g, r, i, and z filters respectively. 
Limiting magnitudes of ~ 26.5, 25.5, 25, and 24 were reached 
in the g, r, i, z filters. 

Several objects are detected within 10 arcsec of the 
QSO, as can be seen in Figure [1] Objects labelled 5, 6 and 
7 are the brightest galaxies in the region, with m;=23.86, 
m,i = 20.90 and rrit — 21.89. Object 6 and 7 appear to be 
interacting. All of these three objects have relatively flat 
spectral energy distributions, typical of actively star form- 
ing galaxies at the redshift of the absorbers. Photometric 
redshifts for these three objects determine z p hot ~ 0.70, 
consistent with these being members of a group. We note 
that the la error on the photometric redshift of object 5 
is ~ 0.11, i.e. the photometric redshift determined from hy- 
perz is only 2a away from the redshift of the higher z system 
at z a b s =0.9281. The extremely blue colors of objects labelled 
three and four in Figure [T] likely exclude these galaxies of 
being at the redshift of either absorber, hyperz determines 
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ID 


A a 


A<5 


m g 


i ii/ 


nt, 


TO- 






ii 


it 











Q0826-2230 



Q1009-0026 
Q1228+1018 

Q1323-0021 
Q1436-0051 



1 






16.18±0.01 


1 T" ?~7 A A Al 

15.74±0.01 


15.50±0 


01 


1 r n c i a ai 

15.3b±0.01 


2" 


1.5 


0.2 


19.1o±0.10 


18.73±0.10 


18.53±0 


10 


18.52±0.10 


3 


—0.8 


6.0 


20.90±0.05 


29.88±0.03 


19.51±0 


03 
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0.2 
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>24.00 


3 
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35 


>24.00 
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-4.5 
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>24.00 
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23.72±0.17 
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10 
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2.8 


-5.5 


22.49±0.10 


21.62±0.06 


20.82±0. 


05 


20.69±0.08 


7 


1.7 


-5.7 


23.24±0.14 
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21.81±0 


08 


21.82±0.14 


8 


5.9 


6.6 


25.43±0.37 
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Table 2. AB magnitudes for the objects labeled in Figure [T] The QSO itself is listed as object 1 in each field. Magnitudes given are 
the isophotal values given from SEXTRACTOR for objects detected above a 3cr significance. Limiting magnitudes are given for objects not 
detected in a given band. a Magnitudes for this object were determined from fitting the profile of the object after subtraction of the 
brighter object. b This object is detected in the PSF subtracted i band image, but the PSF subtraction is not adequate to determine a 
magnitude. 
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Figure 2. Colour combined RGB images of the fields. In all cases, the g, 
Q1228+1018 for which the r, i, and z frames were used. 



and i frames where used in the RGB channels, except for 



photometric redshifts of z p hot = 0.25 ± 0.10, 0.18 ± 0.10 re- 
spectively. 



4 RESULTS 

4.1 Luminosities and K-corrections 

In order to determine the absolute magnitudes of the galax- 
ies observed in these fields, K corrections must be applied 
for an accurate estimate. We have determined i-band K cor- 
rections for these galaxies assuming both a starburst type 
galaxy and an Sa type galaxy, with the templates originating 
from the SDSS galaxy spectra template libraries available at 
http : //www. sdss . org/dr7/algorithms/spectemplates/index 



Although in principle the galaxies hosting the absorbers 
in question could be of early type, the fact that there 



is significant continuum emission at short wavelengths 
apparent from the g band images and flat spectral energy 
distributions for these galaxies indicates that star formation 
is currently ongoing in these systems. Hence, we have 
chosen to include only the templates for these star forming 
galaxies for the K corrections. Assuming an early type 
galaxy spectrum would increase the absolute magnitudes 
of these galaxies. The K corrected absolute magnitudes 
are given in Table |4~T1 assuming that each galaxy is at the 
redshift of the abso rber. The value of M* = —21.59 from 
iBlanton et aD (|2003h was used for determining luminosities. 

In total, the fields of Q1009-0026, Q1228+1018, 
Q1323-0021, and Q1436-0051 all contain L > L* galax- 
hinslwithin 75 kpc projected impact parameters. Could these 



galaxies be unrelated field galaxies that just happen to lie 
near the sightline to the QSOs? This is a persistent question 
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Figure 3. The expected number of galaxies with L>0.1L* in a 
pencil beam of area 100 arcsec 2 as a function of red shift, assuming 
the Sc hechter luminosity function parameters from lBlanton et al.l 
[120031 1. 



when imaging QSO absorbers, which can ultimately only 
be settled with followup spectroscopy to confirm the red- 
shifts of the galaxies. However, the probability of observing 
an L > 0.1L* galaxy in a small projected diameter pen- 
cil beam survey such as this is quite small. Assuming the 
Schechter function p arameters of th e i ban d galaxy luminos- 
ity distribution from lBlanton et all l|2005t ) of M* = -21.59, 
4>* = 0.005 Mpc -3 , and a = —1.00, we have estimated the 
number of galaxies with L > 0.1L* that will be intercepted 
in a random 100 arcsec 2 field, shown in Figure [3] as a func- 
tion of redshift. As can be seen, we expect to intercept less 
than one L > 0.1L* galaxy in a random 100 arcsec 2 field 
for z < 1.5. Hence, the likelihood of each of these fields 
containing bright galaxies not associated with the hosts of 
the sub-DLAs are very sm all. Using th e lumi nosity function 
parameters at z ~ 1 from lllbert et all (|2005h results in the 
same conclusions. 



4.2 Column Density vs. Impact Parameter 

We have compiled the available data from the literature on 
impact parameters for absorption systems that have known 
Nh i at z < 1. These data are given in Table [4] We have 
taken a sub-sample of these data to include only the sys- 
tems that have spectroscopically confirmed galaxy redshifts, 
or secure photometric redshifts. A Kendall's r test was used 
to determine the probability of a correlation between Nh i 
and impact parameter p, resulting in r = —0.23 and a prob- 
ability of no correlation of 0.47. If one includes the data 
points that do not have spectroscopic redshifts of the galax- 
ies, Kendall's r then becomes r = —0.48 with a probability 
of no correlation of 0.05. 

We have preformed a linear regression on these data, us- 
ing survival analysis to include the mixed censored data in 
the fit using Schmitt's Binning method. We find the best fit 
trend line of log N H i =(20.67 ± 0.23) - (0.028 ± 0.014) x p 
when using only the spectroscopically confirmed galaxies, 
and log N H i =(20.86 ± 0.14) - (0.03 ± 0.008) x p when in- 
cluding galaxies without redshift information. 

Figure [4] shows a plot of impact parameter vs. Nh i for 



Figure 4. Log Nh i vs impact parameter for QSO absorbers with 
followup imaging data. Square data points have spectroscopic or 
secure photometric redshift determinations of the galaxies, and 
circular points have no redshift information on the galaxies. The 
shaded region represents the l<r confidence intervals on the best fit 
linear regression line using the Schmitt method to include upper 
and lower limits in the spectroscopically confirmed sample. Also 
overplotted is the 21 cm profile of the LSB galaxy Malinl from 
iPickering et al.Nl994) . and the average 21 c m profile of local star 
formin g galaxies from I S waters et al. | ( l2002h ; | Noordermeer et all 
(|2005h . 



the literature data, and the points from the observations 
in this investigation. For Q1228+1018, and Q1436-0051 we 
give the impact parameter as a lower limit as it is uncertain 
to which galaxy in these fields the absorption is linked to 
without spectroscopic observations of the galaxies. For the 
absorber at z a (, s =0.8426 in Q1009— 0026 we give the impact 
parameter as an upper limit based on a 2" threshold radius 
inside the PSF. Also shown in Figure [4] are the radial 21 
cm profile o f Malinl, a nearby low surface brightness (LSB) 



cm prohle o t Malinl, a near b y low surface brightness (Ls>o ) 
galaxy from IPickering et al.l (Il994|) and t h e aver age Sc type 
21 cm radial profile from Swaters et al.l |2002l ). The data 
in that work are normalized relative to the 21 cm H I ra- 
dius (i.e. R/Rir), so w e have used the median H I radius 
for L~L* galaxies from iNoordermeer et alj (|2005h of 28 kpc 
to scale the profile. The minimal correlation between im- 
pact parameter and Nh i indicates that the QSO absorber 
population stems from a range of environments and morpho- 
logical types, which is also supported in the heterogeneous 
collection of galaxies that comprise the literature sample. 

The random orientations of galaxy disk inclinations in- 
creases the dispersion in the trend between impact parame- 
ter and Nh i- Monte Carlo simulations of p vs Nh i reported 
m (|Chen et all |2005h have a scatter of ~0.5 dex in Nh i 
for a given impact parameter. These simulations were con- 
ducted using idealized smooth H I profiles, using more re- 
alistic simulations that include the patchy and filamentary 
features typically seen in the ISM of the Milky Way and 
nearby galaxies would likely increase the scatter even more. 
A larger sample of DLAs and sub-DLAs with imaging data 
at z < 1 would allow for a direct comparison of objects in 
terms of their impact parameters an luminosities. 
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Table 4. a n = no redshift information, p = secure photometric 

redshift, s = spectros copic redshif t 

REFE REN CES: (1) IChen et all J2003T). (21 iBorthakur et al.l 
d2010ll. (3) iTurnshek et alj j200ll ) , (4) jChen et all (120051). (5) 



lEllison. Kewlev fe Mallen-Ornelaaf <2005h . (6 
1120051). (7) lLanzetta et al.l h997l ), (8) 



.Christensen et all 
Rao et al .1 <2003h. 



9) iKacprzak. Murphy fc ChurchiliT j2010h.(10) lie Brun et all 



(13) IBurbidge et al.l U99a), (1 4) lUhun et alj (I2U11I ) , (15 ) 
Work, (16) iLacv et all l|2003l ). (17) IPeroux et all (201Cl) . 
iHewett fc Wild I ( 120071) 



(11) ISteidel 



et all Jl994l ) (12) iGharanfoli et all l|2007h 



jl996l), (1 4) jchun et all ||2010| ) (15 ) This 



(18) 



5 DISCUSSION AND CONCLUSIONS 

In our previous spectroscopic observations of sub-DLAs 
at 2 < 1.5, we have shown that this population is more 
metal rich than their DLA counterparts. With the ob- 
se rved correlation betwe en luminosity and metallicity seen 
in iTremonti et al. I (|2004h . one would suspect that the host 
galaxies of the metal rich DLAs in our sample should be 
L>L*. We have preformed deep imaging of the fields of five 
QSOs that harbour metal rich sub-DLAs in their spectrum 
with the SOAR telescope and the SOI camera. In four of 
the five fields, luminous galaxies are observed, and one QSO 
(Q0826-2230) is likely gravitationally lensed by the interven- 
ing absorbing galaxy. 

Is it possible that the candidate galaxies are directly 
atop the QSOs? The PSF subtraction of Q1323-0021 suc- 
cessfully reveals the galaxy detected in I Chun et al I2010I . and 
no other bright galaxies are detected in the PSF subtracted 



images so we feel that this is unlikely. There are several rea- 
sons to suspect that the galaxies at large impact parameters 
are in fact the host galaxies of these sub-DLAs. The ever pos- 
sible scenario of a LSB galaxy at low impact pa rameter is 
unlikely, as the low metallici ties of LSB galaxies (|McGaughl 
Il994l ; iHaberzettl et al1l2007l ) exclude them as possible can- 
didates for these high metallicity sub-DLAs. The probability 
of observing any interloping field galaxy with L>L* in a pen- 
cil beam survey such as this is small, as is shown in Figure 
[3] and the photometric redshifts for the galaxies are consis- 
tent with these galaxies being at substantial redshifts. Aside 
from the case of Q1323— 0021 in which the impact param- 
eter is ~9 kpc, in the fields of Q1009-0026, Q1228+1018, 
and Q1436— 0051 the impact parameters of all the galaxies 
in the fields is >35 kpc. 

These large impact parameters in several cases cer- 
tainly raise the question as to the origins of the absorb- 
ing material. Could the material be ordinary interstellar 
gas residing in the extended disk of the host galaxies? The 
H I envelope of gas rich galaxies typically does extend to 
much greater dista nces than the optical images suggest 
l|Walter et all 120081 ) , although not often to great enough 
radii to produce sub-DLA H I column densities at these 
impact parameters. Above a certain critical surface den- 
sity, l ocal pockets of gas are unstab le to gravitational col- 
lapse (|Kennicutdll98alToomre!ll964l ). The low H I column 
densities of these sub-DLAs are well below the typical H 
I colu mn density threshold of ~ 10 20 ' 5 cm -2 of iKennicuttl 
(|l989f ). hence star formation is unlikely to occur in these 
regions. Observations have shown some amount of star for- 
mation at large galactic radii in nearby galaxies, however 
the star formation complexes are sparse with a low fill- 
ing factor, making a r andom sightline passing through such 
regions improbable ((Ferguson. Gallagher fc Wvsei Il998 



iGoddard. Kennicutt fc Weberl I2010T ). Merging galaxies do 
however often show high levels of star formation at la rge 
galactic radii ([Larson fc Tinslevlll978l ; ISmith et alj|2010h . 

Similarly, the abundances as derived from the nebu- 
lar lines in the outer regions (R>R2s) of nearby late type 
galaxies are typically ~ 0.1 Zq even where the abundances 
are roughly solar in the inner regions of the galaxy. These 
metallicities are similar to that of many DLA systems, but 
well bel ow the abundances we see in the sub-DLAs in this 
samp le (|Ferguson. Gallagher fc Wvsei Il998l ; iBresolin et all 
l2009h - The impact parameters we see in these absorbers is 
several times greater than R25 ~ 10 — 15 kpc typical of late 
type galaxies. 

The relatively flat spectral energy distributions for the 
galaxies in the fields of Q1009-0026 and Q1436-0051, for 
which we have four filter imaging are indicative of star 
forming galaxies. In the case of Q1436— 0051, there is clear 
indications of interacting galaxies which are typically ac- 
tively forming stars. We feel that a more likely explana- 
tion for the origins of this gas in light of these imaging 
observatio ns is that it is the signature of outflowing galac- 
tic w inds (jBouche et al.|[2007l ; lKobavashi. Springel fc White! 
l2007h . The high metallicities of these systems along with the 
large impact parameters suggest that this gas is the result 
of recent bursts of star formation in the galaxies or possibly 
tidally stripped material in the interacting galaxies. 

The mean Augoof sub-DLAs as a population is ~125 
km s _1 , slightly higher than the mean for DLAs of 
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103 km s" 1 (jKulkarni et alj l20ld L The sub-DLAs in 
this work have a slightly lower mean velocity width than 
the total sub-DLA population of Avgo=107 km s _1 . 
The complex velocity profiles seen in some absorption 
line systems have of ten been attributed to merg i ng or 
interacting galaxies jPetitiean. Srianand fc Ledoux |2002| ; 



iQuast. Reimers fc BaadeT 20081 ) . The profiles for both the 
systems in the spectrum of Q1436— 0051 are not unusu- 
ally broad, each extending ~100 km s _1 . If the interacting 
galaxies seen to the south of this field are indeed responsi- 
ble for the absorption line system, then at least in this case 
interacting galaxies surprisingly do not seem to produce a 
complex velocity structure. At smaller impact parameters 
the profiles may show signs of broader velocity profiles. The 
inclination of the disk to the QSO line of sight possibly plays 
a large part in the kinematic width of the absorption profiles. 
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It is i nteresting to note that alth ough the recent com- 
pilation of iNoterdaeme et al.l (|2008rJ ) focused on DLA sys- 
tems (77 total systems were included in the sample, with 
only 7 sub-DLAs and all had log Nh i >20.0), one of the 
systems with a high molecular fraction was in fact a sub- 
DLA. Similarly high molecu le rich sub-DLAs have also been 
seen in other investigatio ns (|Quast. Reimers fc Baadell2008l ; 
INoterdaeme et aT, 2008a). It would be interesting to expand 
the number of H2 measurements in sub-DLAs to see if they 
are perhaps also have high molecular fractions, and thus fa- 
voring the scenario where the lower H I column densities are 
due to the conversion of gas from neutral to molecular. 

Although the number of DLAs or sub-DLAs that have 
been observed with high quality spectra has increased, there 
is still a very small number of galaxies that have been con- 
firmed at the redshift of the absorber with followup spec- 
troscopic measurements. Even fewer DLAs have spectra of 
the host galaxy of high enough S/N that can be used to 
determine the properties of the galaxy's stellar population 
such as the star formation rate (SFR), masses and star for- 
mation histories via spectral template fitting and emission 
line diagnostics. 

Followup spectroscopy is clearly necessary to link the 
properties we see in absorption such as the metallicity, kine- 
matics, and abundance patterns with the properties of the 
stellar populations of the host galaxies. DLAs and sub-DLAs 
are likely to arise in a range of environments, a larger sam- 
ple with imaging and spectroscopy of the host galaxies would 
reveal any potential dichotomy in the populations. The sub- 
DLA systems in this sample seem to arise in a range of 
environments themselves, from the inner region of a possi- 
ble early type galaxy in the case of the sub-DLA in SDSS 
J1323— 0021, from the periphery of a luminous star forming 
galaxy in SDSS J1009— 0026 and from the outskirts of an 
interacting pair of galaxies in SDSS J1436— 0051. However, 
in all cases we observe luminous galaxies that are the likely 
host galaxies of these metal rich absorbers. A large sample of 
both spectroscopic and photometric measurements of DLA 
and sub-DLA host galaxies will allow for studying trends be- 
tween luminosity, impact parameter, and SFR. Higher spa- 
tial resolution imaging with the HST would also allow for 
detailed information on the morphology of the absorbing 
galaxies. 
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